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ABSTRACT 


Integrators are included in the feedback loop of a control system to eliminate the 
steady state errors in the commanded variables. The integrator windup problem arises if 
the control actuators encounter operational limits before the steady state errors are driven 
to zero by the integrator. The typical effects of windup are large system oscillations, high 
steady state error, and a delayed system response following the windup. In this study, 
methods to prevent the integrator windup are examined to provide Integrator Windup 
Protection (IWP) for an engine controller of a Short Take-Off and Vertical Landing 
(STOVL) aircraft. An unified performance index is defined to optimize the performance 
of the Conventional Anti-Windup (CAW) and the Modified Anti-Windup (MAW) 
methods. A modified Genetic Algorithm search procedure with stochastic parameter 
encoding is implemented to obtain the optimal parameters of the CAW scheme. The 
advantages and drawbacks of the CAW and MAW techniques are discussed and 
recommendations are made for the choice of the IWP scheme, given some characteristics 


of the system. 



Chapter 1. 


Introduction 

Flight control systems for modem tactical/fighter aircraft are to be designed for 
enhanced flight maneuvers such as Short Take-Off and Vertical Landing (STOVL) and 
high angle of attack flights. The forces and moments generated by the conventional 
control surfaces are usually insufficient for successful completion of such maneuvers 
under the flight conditions considered, and hence are augmented with the forces from the 
propulsion system. The engine and its subsystems also play an important role for aircraft 
with tilt-rotors and for aircraft with engine steering control. Thus the propulsion system 
of a modem aircraft is not being used to just overcome the atmospheric drag, but also to 
control the airplane. Hence a need arises to control the engine variables as precisely as 
possible. 

1.1 Integrator Windup Problem 

One of the difficulties in controlling the engine variables arises due to the 
integrator windup, caused due to the control actuators encountering the operational limits 
when integrators are used in the feedback loop of the control system. Operational limits 


i 
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include both physical limits on the actuator motion as well as the user imposed limits on 
certain variables of the system to avoid extreme operational conditions. For instance, 
limiting the maximum fuel flow to the engine to avoid high temperature damage of 
turbine blades is an operational limit on the fuel flow. Integrators are essential in the 
feedback loop to drive the steady state errors in the commanded variables to zero. If a 
control actuator encounters the operational limit before the steady state errors reach zero, 
the errors cannot be driven to zero by the controller. The integrator would continue to 
integrate this non-zero error and hence the integrator’s output would build-up to a very 
high value. This phenomenon is referred to as integrator windup, and results in increased 
system oscillations, degraded system performance, and in some cases even an unstable 
system [1], The windup effects also arise if the dynamics of the controller is relatively 
slower them that of the error when the operational limits of the system are encountered. 
Following this period of windup, the controller’s response to new command inputs might 
be poor because the integrators must first unwind prior to attempting to drive the steady 
state errors, due to new command inputs, to zero [2], 

In the early phases of feedback control, practicing engineers had developed a good 
intuition on the effects of integrator windup and methods for handling it. Often, the 
integral action was achieved due to the motion of the control actuator itself [3]. In such 
cases, integration stops when the actuator encounters the physical limit, and hence the 
windup is avoided. These methods had little theoretic foundation but worked well for 
Single Input-Single Output (SISO) systems that are simple to analyze. However, modem 
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day systems are mostly multivariable, and such Multiple Input-Multiple Output (MIMO) 
systems are quite cumbersome to analyze due to the complexity of the system. In a 
MIMO system, the windup problem is further aggravated if there exists a strong coupling 
between the control variables of the system, which degrades the system performance even 
if a single control actuator encounters the operational limits. Also, if the MIMO system is 
sensitive to direction changes of control vector, then the performance of the system is 
further deteriorated [4], While much progress has been made in the development of 
multivariable control theory, some important issues like the integrator windup problems 
have not been addressed yet in a systematic manner. 

It might seem to be a rational approach to handle the windup problem at the 
design stage of the controller by including the plant limits [4]. However, this approach is 
quite tedious and the resulting control law is often very complicated. Also, the 
nonlinearities of the actuator are not always known apriori. Hence, the design of the 
controller is often based on linear theory by neglecting the limits. An extra feedback 
compensation is then added at the control implementation stage to take the practical 
limits into consideration. As this compensation aims to diminish the effects of windup, it 
is referred to as anti-windup or Integrator Windup Protection [5]. This additional 
compensation leaves the original linear behavior unchanged but provides graceful 
degradation of system performance when actuator limits are encountered. 
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1.2 Approaches to providing IWP 

Most of the IWP schemes modify the error between the commanded input and the 
plant output to achieve the anti-windup action. Typically, this error is reduced by adding 
a non-zero factor that is proportional to the extent of the windup. This decreases the 
magnitude of the error that would be integrated by the integrator and hence diminishes 
the effects of windup. This principle is also called back-calculation since the controller 
states are back-calculated such that the output of the controller is at the actuator limit [3], 

In the Conventional Anti-Windup (CAW) scheme, the difference between the 
actuator limit and controller output is multiplied by a set of gains (referred to as IWP 
gains) and is added to the above error term when the actuators are limited. The CAW 
method works well for a SISO system, but for a MIMO system, anti-windup provided by 
the CAW scheme may not be adequate [6]. Also, the CAW scheme does not maintain the 
direction of the controller’s output vector, and hence is not suitable for systems that are 
sensitive to changes in control vector direction. 

In the Modified Anti-Windup (MAW) approach, controller states are multiplied 
by a windup factor before being added to the above described error term. The windup 
factor is such that the direction of the controller output vector is maintained [1 ] and hence 
this scheme is suitable for plants that are sensitive to direction changes of the control 
vector. The main drawback of the MAW scheme is that it has few design variables and 
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hence we do not have sufficient degrees of freedom to design the IWP for time varying 
actuator limits. Consequently, the performance of the MAW scheme suffers when the 
system encounters a limit different from the one for which it was designed. 


1.3 Optimization Methods 

The choice of the optimization method used to determine the set of optimal 
parameters is an important step in the overall solution process. Genetic Algorithms (GA) 
have been successfully used in the past to solve a number of non-linear optimization 
problems [7], As the GAs do not require any gradient information to aid their search, they 
can be used even for cases where the derivatives are difficult to obtain. However, GAs 
require information about the range in which parameters of the optimization problem are 
expected to lie. This handicaps a regular GA approach in cases where we do not know the 
expected range of parameters apriori. Also, for problems with large number of 
optimization parameters, the regular GA needs a large population size to find the optimal 
solution, and this calls for enormous computational resources [8]. In this study, an 
improved version of the stochastic GA proposed in Reference [9] is implemented to solve 
the IWP optimization problem. This modified GA with a stochastic encoding structure 
overcomes the shortcomings of the regular GA in that it does not need the exact 
parameter ranges, and it converges to the global optimum quickly with a small 
population. This search technique starts with a given search region and as the GA 
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population evolves, the search region is modified to capture the global optimum and 
finally to converge to it. The accuracy of the GA solution is verified by carrying out the 
optimization using the optimization module supplied with MatrixX [10]. 

1.4 Study Objectives 

In this study, we develop an unified performance index to provide anti-windup 
compensation using the CAW scheme. This performance index is optimized to provide 
IWP for the engine control system of a STOVL aircraft. The IWP is implemented on the 
linear model of the engine control system developed at NASA Lewis Research Center. 
The objectives of this study include : 

• Studying the CAW and MAW schemes for providing IWP and implementing 
them to example problems to show their merits and demerits. 

• Development of an unified performance index for the design of IWP. 

• Development of an improved stochastic Genetic Algorithm (GA) method to 
obtain the optimal IWP gains. 

• Application of the CAW and MAW schemes to provide IWP for the engine 
controller of a STOVL aircraft. 

• Analysis of the results obtained using the CAW and MAW schemes. 

• Conclusions regarding the suitability of CAW or MAW for a given control 


system. 
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1.5 Thesis Organization 

This thesis is organized as follows: In chapter 2, the previous literature on the 
IWP methods is detailed followed by the optimal parameterized IWP formulation for a 
generic control system. In chapter 3, we examine the optimization techniques used to 
obtain the optimal IWP gains. In chapter 4, we implement the IWP schemes for the 
engine control system of a STOVL aircraft, and discuss the results obtained. The 
conclusions are presented in chapter 5, followed by the recommendations for the choice 
of the IWP scheme, given some system characteristics. 



Chapter 2 


Optimal Parameterized IWP Formulation 

Most natural and man made systems are inherently non-linear, though to a varying 
degree. It is well known from the non-linear control theory that analysis and design of 
non-linear control systems is quite complex and mathematically unwieldy. Hence, for the 
ease of analysis we linearize the non-linear system about a particular operating point and 
the resulting linear system provides an acceptable solution as long as we stay close to the 
operating point However, in most practical cases, large disturbances acting on the system 
and changes in operating points cannot be avoided, thereby forcing us to face the non- 
linearities temporarily. Examples of such non-linearities include actuator limits, 
parameter variations due to changes in operating point of the system, backlash in gears or 
valves [11], hysteresis [12] etc. 

The effect of feedback control is to reduce the system sensitivity to such external 
disturbances, to minimize the effect of plant parameter variations, and to modify the 
system dynamics to a desired form. If in addition, the steady state errors on the 
commanded variables are to be driven to zero, it can be shown that the feedback loop 
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should have an integral control law. However, when integrators are used in the feedback 
loop, limits on the actuator can cause the integrators to windup thereby deteriorating the 
performance of the system. Further, the response of the system following the period of 
windup could be sluggish until the integrator unwinds. Integrator Windup Protection 
(IWP) is hence included in the feedback loop to minimize the effects due to windup. In 
this chapter, first we briefly review the available literature on IWP methods. The IWP 
methods we consider for further study, the CAW and MAW schemes, are then described 
in detail followed by the description of the performance index that conforms with the 
requirements of IWP. We then apply these IWP schemes to two example problems to 
illustrate the effect of directional sensitivity of the system on the choice of the IWP 
scheme. 

2.1 Literature Overview 

Early efforts to handle integrator windup were mainly focused on SISO systems. 
According to the Back-Calculation principle given in Reference [3], when the controller 
output exceeds the actuator limits, the integral is re-computed such that its new value 
gives an output at the limit. 

Reference [1] examines two methods based on back-calculation principle to 
provide windup protection for a turbofan engine control system. The first method is an 
extended version of the Conventional Anti-Windup (CAW) scheme for multivariable 
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controllers. In the Conventional Anti-Windup method, the actuator error is fed back 
through the IWP gains as an additional contribution to the state derivative calculations. 
The MAW scheme uses a scalar windup factor to modify the magnitude of the control 
vector while maintaining its direction. The calculation of this windup factor is a non- 
linear function of the controller outputs and the actuator limits. The CAW and MAW 
schemes are described in detail in a later section. 

Reference [6] points out that the directional sensitivity of a MIMO system plays a 
key role in the design of IWP. Examples are given in reference [6] to show the that the 
CAW scheme does not provide satisfactory protection against windup for directionally 
sensitive systems. The structured singular value (p) analysis is used to determine the 
directional sensitivity of the CAW system. We include these examples at the end of this 
chapter to illustrate the effect of directional sensitivity of the system on the performance 
of the CAW system. 

Reference [4] includes the issue of directional sensitivity in the IWP design 
requirements. The structured singular value (p) is used in [4] to determine the sensitivity 
of the system to direction changes of the control vector. We adopt this measure later in 
our performance index definition (section 2.6) to minimize the directional sensitivity of 
the CAW system. It is pointed out in reference [4] that, if other IWP schemes do not 
satisfy the p requirements, the MAW scheme could be used to provide anti-windup 


compensation. 


Reference [13] implements an intelligent limiter based on back-calculation 
principle for a second order plant and carries out the stability analysis of the limited 
system using the describing function method and Nyquist stability theorem. Reference 
[11] implements a modified back-calculation principle for windup protection for cascade 
controllers, where a limit at secondary actuator may cause windup in the primary 
controller. 

In Incremental Algorithms approach [3], the rate of change of the control signal is 
first computed and then fed to an integrator. Integration is stopped whenever the output 
exceeds the actuator's operational limits. Reference [5] implements the discrete-time 
implementation of Incremental Algorithm for a PID controller. 

Conditional Integration technique is implemented in [5] to avoid windup of a PID 
controller. According to this technique, when the controller works in the linear region, i.e. 
when the controller output is not limited, the error term is integrated by the integrator and 
when the actuator is limited, the integration is stopped. Reference [3] introduces the 
notion of a proportional band and integrates the error only if the predicted process output 
is in the proportional band. The proportional band is the interval in which the process 
output lies when the control variable is varied within the control limits. Reference [14] 
implements an anti-windup method based on conditional integration to provide windup 
protection for a digital multivariable controller. It, however focuses on a particular system 
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with two inputs and two outputs and the implementation is difficult to generalize for 
systems with more inputs or outputs. 

In the conditioning technique approach proposed by [12], a realizable reference 
signal is input to the system instead of the actual reference input. The realizable reference 
signal is computed such that if it is applied to the controller, the controller output would 
have been at the actuator limit. Reference [4] shows that this technique might provide 
poor anti-windup performance for some cases. 

In the observer based approach, the controller states have the physical 
interpretation as the estimates of plant states [4], Hence the objective here is to design 
such that the controller states assume the correct estimates of the plant states regardless of 
plant limitations. If the controller has full access to the plant states, this can be achieved 
by a simple state feedback. In the case where the full plant states are not accessible, an 
observer is constructed to supply plant states which tire used to provide the feedback. 

Internal model control (IMC) based anti-windup technique requires an exact plant 
model for closed loop stability when the actuator encounters operational limits. Reference 
[4] shows that there are no inherent properties of IMC that provide robustness to diagonal 


input uncertainties. 
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Most of the anti-windup methods described above are applicable only for SISO 
systems. Some of the anti-windup examples have been applied to MIMO system 
examples, but they are highly dependent on the problem being solved. They also lead to a 
complex setup when the number of inputs/outputs increase. Hence a need is felt to derive 
a simple anti-windup solution that is applicable to a generic MIMO system regardless of 
the number of system inputs/outputs. We begin by stating the general requirements of an 
IWP scheme. 


2.2 IWP Design Requirements 

Integrator Windup Protection (IWP), when included in the feedback loop should 

prevent the integrator’s output from building up when the actuator encounters a limit 

Some of the important design requirements for an IWP scheme are[l ]: 

1 . A limited actuator must be observable. 

2, IWP should be memoryless and should not contribute to the control system when the 
limits are not encountered. 

3 IWP should provide closed loop stability for all possible actuator limit combinations 
within the system’s operating envelope. 

4 IWP should attempt to maintain system performance for all possible actuator limit 
combinations within the system’s operating envelope. If the system performance 
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cannot be maintained, IWP scheme should provide a smooth, stable transition to a 
minimally degraded operating point. 

5. IWP should provide smooth transfers between the unlimited and limited actuators, 
while providing accurate tracking of the limited actuators. 

6. IWP should minimize the sensitivity of system to direction changes of the control 
vector due to actuator limits. 

We now implement an IWP scheme for a generic control system, which, together with the 
structure of IWP and the definition of performance index would fulfill these 
requirements. 

2.3 A General IWP Structure 

A generic structure of a control system with actuator limits and a IWP structure is 
shown in Figure 2.1. In this figure. G(s) is the plant, the ''Limits' block represents the 
limits on the actuator, K(s) is the controller and IWP is the anti-windup compensation 
added to the controller as shown. The reference command to the control system is the 
controller output is ‘« c ’, the plant inputs are limited by the limits ‘ u L ' and the plant 
outputs are ‘z\ The error V between the reference inputs and the plant outputs is fed to 
the controller K.(s). Also, the information about the plant inputs, l u L c ' , is supplied to the 
controller in order to implement the IWP. It is necessary that the IWP scheme included in 
the controller structure stabilizes the controller and the closed loop system. If there is no 
error between the plant inputs and the controller outputs (i.e. when the limits not 
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encountered), the IWP gains should be ineffective. When the plant inputs are different 
from the controller outputs (as in the case of actuator limitation), IWP should reduce the 
magnitude of error between the plant outputs and the reference inputs so that the windup 
effects are diminished. We assume that for all systems considered in this study, the 
actuator is observable. This satisfies the requirement # 1 of the IWP design requirements 
stated in section 2.2. 


Controller 


Limits 


Plant 



Figure 2.1 General IWP Structure 


We now describe two IWP methods that satisfy the requirements stated in section 2.2, 
and are applicable to a generic MIMO system. These IWP methods, namely the 
Conventional Anti-Windup method and the Modified Windup method, are based on the 


back-calculation principle. 
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2.4 Conventional Anti-Windup (CAW) Scheme 

In the CAW method, the controller states are back-calculated by using actuator 
error feedback such that output of the controller is within the actuator limits [1). The 
actuator error is the quantity by which the output of the controller exceeds the actuator 
limits. An implementation of the CAW scheme is shown in the Figure 2.2. As seen from 
the figure, if the controller output is within the actuator limits, then the actuator error 'e u ' 
is zero and the CAW gains are ineffective. Thus the CAW scheme is memoryless and 
hence satisfies the second IWP design requirement. When the controller output exceeds 
the limits, the actuator error i e u ' is fed back through the gain matrix ‘A’, so that the 
controller output remains within the limits. The terms in the gain matrix ‘A’, referred to 
as IWP gains, dictate the stability and performance of this system protected against 
windup. These are obtained such that they satisfy the stability requirements and optimize 
the system performance. 

The equations for the CAW implementation can be given as[l] : 


e 

y_ 


u = C x + D 
c c c c 


e 

y 


*c 


= A x 
c c 


+ B 


Original Dynamic Controller : 
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Controller with IWP : 


= Vc + B c 


+ A le. 


= (a c -U.C c )x c+ [b c -M.D c \ 


+ A Lu 


L 


In the above equations, [A c , B c , C c , D c ] are the system matrices for the controller, x c is the 
state vector of the controller, [e y] T is the controller input, u c is the output vector of the 
controller, u c L is the vector of actuator limits, e u is the actuator error vector, A is the 
constant IWP gain matrix and L is a diagonal matrix that represents the actuator being 



Figure 2.2 CAW Scheme Implementation 

limited. If L=0, the nominal control system is obtained. For a two actuator system, 
L=diag(l,0) represents first actuator being limited, L=diag(0,l) denotes that the second 
actuator is limited and L=diag(l,l) denotes that both actuators have encountered the 
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limits. As seen in the above equations, the stability of the CAW controller can be 
obtained from the matrix ( A c -ALC c ) for various combinations of L and the chosen IWP 
gains A. 

When a control variable encounters the actuator limit, the CAW scheme truncates 
the variable such that the control is within the limit. This is illustrated in Figure 2.3 [1], in 
which a two actuator control system is considered. In this figure, u cl and u c2 are the two 
control variables, u u and u L2 are the actuator limits for these control variables, and d, is 
the direction of the unlimited control vector ( u c ). As seen from the figure, the control 
vector modified by CAW scheme Uc AW has a direction d 2 different from the unlimited 
control vector. Some systems are sensitive to changes in the direction of control vector 
and for these systems, the CAW scheme could lead to a poor closed-loop performance. 
Hence, for such systems, we must ensure that the CAW gains minimize the sensitivity to 
the control vector direction changes. A distinct advantage of this scheme is that the 
number of design variables (terms in the gain matrix A) is equal to the product of number 
of actuators and the number of states of the controller. Hence we have sufficient degrees 
of freedom to design the IWP for an acceptable performance level. 

2.5 Modified Anti-Windup (MAW) Scheme 

In this anti-windup method, modified control vector is obtained by scaling down the 
unlimited control vector such that the control variables are within the actuator limits. 
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When the control vector is scaled down by a scalar factor, the direction of the control 
vector is maintained. In Figure 2.3, is the control vector modified by the MAW 
scheme and as seen from figure, it has the same direction as the unlimited control vector. 
Since this scheme maintains the direction of the control vector, it is suitable for systems 
sensitive to direction changes of control vector. Figure 2.4 shows the implementation of 
the MAW scheme. In this figure, a is the scalar windup factor that denotes the extent of 



i 


Figure 2.3 Control Vector Directions for CAW and MAW schemes 


windup and p is a design variable. The windup factor a is defined as 
a (t) = l if u c is within limits, otherwise 


a(/) = min 

i 
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where u L c is the limited control vector (plant inputs), u c is the unlimited control 
vector (controller outputs), and subscript / denoting the / component of the vector The 
design variable P is to be determined such that the closed loop system is stable and the 
performance of the IWP system is close to that of the nominal system. As seen from 
Figure 2.4, the nominal control system is active as long as a = 1. When a control actuator 
encounters the limit, a*l, and the nominal controller is modified by the two additional 
blocks shown in the figure. The state vector x c is fed back after multiplication by the 
factor P(a -l) so that the output of the controller is eventually brought within the actuator 
limits. 
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Figure 2.4 MAW Scheme Implementation 
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The equations that implement the MAW scheme can be given as [11: 


*c=K+P(a-l)l)* c+ B c 


u = C x + a D 
c c c c 


y 


The system matrix (A + P(ct - 1)/) in the above equation can be rewritten using 
similarity transformation as 

r'(A + (3 (a - l)/)r= r'AT+(}(a - 1)/ 


Since 0< a <1 and p > 0, it follows that MAW scheme always shifts the real part of the 
controller eigen values to the left, thereby making the controller more stable. 

As described earlier, MAW scheme works well for systems sensitive to control 
vector direction changes. This is illustrated by an example provided at the end of this 
chapter. However, since the design parameter of this approach is a single scalar variable 
(P), we do not have a wide degree of freedom to tune the system for an optimal 
performance. Also, the optimal performance might require a high value for the gain P that 
is not be feasible to implement in practice. The non-linearity associated with the MAW 
scheme while computing the scalar windup factor might degrade the system performance 
though the effects of windup are eliminated. 
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2.6 .Performance Measures for IWP 

In this section we discuss the performance index chosen to optimize the 
performance of an IWP scheme. The design parameters of the IWP scheme are to be 
determined before we could use it to provide windup protection. The parameters have to 
be chosen such that the stability of the closed loop IWP system is maintained and the 
IWP system performance is close to the nominal system performance. The performance 
requirements can be achieved by defining a suitable performance index for the IWP 
system such that when this performance index is optimized, the performance of the IWP 
system is close to that of the nominal system. In Reference [2], the root mean square 
(rms) error between the nominal and IWP system due to a zero mean, unit variance white 
noise input signal, is considered to be the performance measure for the IWP scheme. The 
IWP gains are parameters of the optimization problem, and hence can be determined 
using an appropriate optimization technique. Stability requirements are imposed as 
constraints while optimizing the performance of the IWP system. If a set of parameters 
result in an unstable system, a poor performance index is assigned for that set of gains 
and we proceed to design the IWP with a different set of parameters. Hence we ensure 
that the IWP gains result in a stable controller and a stable closed loop system before 
attempting to evaluate the system performance. This is to conform with the IWP design 


requirement U 3. 
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The IWP optimization structure used to obtain the performance index in [2] is 
illustrated in Figure 2.5. As shown in the figure, the nominal system and the IWP system 
are augmented with signal conditioning blocks that provide appropriate scaling and 
frequency weighting for the inputs and the errors. We adopt the performance definition 
structure for our study from reference [2]. The data for the signal conditioning and 
frequency weighting blocks were obtained from the NASA Lewis Research Center. The 
explanation for each of the blocks is given below: 

1 The “Command Loop Shaping” block consists of a loop scale factor and a first order 
lag for each of the controller loops. The external commands Z c and U L are white noise 
signals that are scaled and filtered by the command loop shaping block such that the 
inputs to the nominal and IWP systems have appropriate magnitude and frequency 
spectra. 

2. The “Performance Error Weighting” block weights the performance errors between 
the nominal and IWP systems such that the low frequency (steady state) errors are 
given a higher weighting. This choice of weighting can be accomplished by choosing 
a weighting function that has a large magnitude at low frequencies and drops off to 
small magnitudes at high frequencies. This block also has a scale factor that is just the 
inverse of the scale factors used in the “Command Loop Shaping” block. 

3 The “Actuator Position Weighting” block weights the actuator position errors so that 
the limited actuator is tracked accurately. This ensures smooth transfers between 
limited and unlimited actuator. This satisfies the IWP design requirement # 5. 
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In Reference [2], the performance index is solved as a white noise covariance 
optimization problem. According to this idea, the rms value of the performance errors and 
actuator position errors due to zero mean, unit variance white noise signal inputs is 
minimized. It can be shown that this rms norm is equivalent to solving the LQG problem 
with the performance index given by : 

J 2 = E{ lim - j(W z Z e 2 + Wy U 2 e )dt} 

r-»oo / Q 

where the weights W z and are chosen to penalize the performance errors and the 
actuator position errors respectively. This choice of the performance index ensures that 
IWP system performance matches closely with the nominal system performance, and 
ensures accurate actuator tracking. 



Figure 2.5 IWP Optimization Structure 
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The system pv-Tormance for all actuator limit combinations is ensured by carrying 
out the optimization for all possible limit combinations. Thus, the optimization structure 
and the rms performance measure adopted from reference [2] satisfy the IWP design 
requirement # 4. 

However, when the CAW gains are obtained by optimizing the rms performance 
measure as in reference [2], the closed loop IWP system might be directionally sensitive. 
In that case, the CAW scheme results in a deteriorated performance when actuator limits 
are encountered. Hence, we have to include the directional sensitivity measure in the 
performance index such that the directional sensitivity of the CAW system is minimized. 

The directional sensitivity is the sensitivity of the system to the direction changes 
of the control vector, caused when the actuators encounter the operational limits. The 
directional sensitivity of the IWP system can be obtained by treating the actuator limits as 
uncertainties (with uncertain limits, the direction of control vector is uncertain). We can 
then reformulate the IWP system such that the limits on the system appear as diagonal 
input uncertainties to the system [6] as shown in Figure 2.6. Now the problem is reduced 
to determining the sensitivity of the system to structured uncertainties (in this case, the 
uncertainties have a diagonal structure) which is addressed by the structured singular 
value (p) theory. Thus robustness of the IWP system to changes in control vector 
direction is given by its structured singular value (p). 
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Figure 2.6 Limits as Diagonal Input Uncertainties 


Structured singular value [15] can be defined as the inverse of the size of the 
smallest perturbation. A, that destabilizes the closed loop system shown in Figure 2.6 
Alternatively, if we are given the perturbation A, we can determine whether the closed 
loop system is robustly stable for these perturbations from the structured singular value of 
the system. If the IWP system can be represented as G(s) with the partition structure as 
shown, then the closed loop system of Figure 2.6 is robustly stable for all perturbations 
within the specified bounds if and only if [15] : 

sup|i[G 22 (y<o)] <1. 

0) 

This states that the closed loop IWP system is robustly stable only if the 
maximum structured singular value of the transfer function between e u and U L C ,over all 
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frequencies is less than unity. Here, U L C is the controller commanded control, and e u is the 
quantity by which the controller commanded control exceeds the actuator limits. In this 
analysis, the actuator limits form the perturbation matrix A, as shown in the figure. This 
perturbation structure is diagonal and a 10% magnitude bound on these perturbations is 
assumed. 

While including the directionality in the performance index for the CAW system, 
we first obtain the CAW gains that minimize the rms norm of the optimization structure 
presented earlier. Once we know the minimum rms norm (J 2m j n ), we can minimize the 
structured singular value of the system while imposing a constraint on the rms norm, so 
that the resulting CAW gains minimize the directional sensitivity of the system, and 
provide a rms norm close to J 2m j n - This can be equivalently stated as : 

min { supp[G 22 0'«)] } 

0) 

subject to J 2 < kJ 2mm 

where we minimize the maximum structured singular value of the system transfer 
function between e u and U L C subject to a rms norm constraint. The factor ’k’ is a scale 
factor that denotes the performance degradation we are willing to tolerate in order to 
obtain a system less sensitive to the control direction. Hence, with this performance index 
definition, we satisfy both the performance requirements (requirement # 4), and the 
directional sensitivity requirement (requirement # 6) stated in section 2.2. It should be 
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noted that the MAW scheme maintains the control vector direction, and hence we need 
not include the ji calculations for the MAW scheme. Thus the 1WP schemes discussed 
here (CAW and MAW), together with the performance index definition satisfy all the 
general IWP design requirement stated in section 2.2. 

We conclude this chapter by implementing the IWP schemes on two example 
problems presented in [6]. These examples illustrate the role played by the directional 
sensitivity of the limited system on the choice of IWP scheme. 


Example 1: 

The first example we study is a nominally stable 2 input 2 output system that is 
represented as : 

System A : 


4(0.1 + 5 ) , 

P = — -/T 1 

s 


with 


4 5> 

3 4 J’ 


K = 


1 


4(0.1 + s) 


R 


In the above representation, P and K represent the plant and the controller 
respectively. The schematic of the closed loop limited system is presented in Figure 2 7 
In this figure, d is the desired input to the plant, y is the plant output and e is the error. 
The limits on the actuator are assumed to be [-1,1] as shown in the figure. The response 
of the nominal and limited systems for a step input of amplitude [0.61 0.79] T is presented 
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in Figure 2.8. The overshoot and large oscillations exhibited by the limited system are the 
effects of windup. Here the windup effects are caused because the dynamics of the 
controller is relatively slower than that of the error when the operational limits are 
encountered. 

Figure 2.9 shows the implementation of CAW scheme for this system. As seen in 
this figure, the actuator error is fed back to the controller through a set of CAW gains. 
Since we are considering the examples to illustrate the effect of directional sensitivity, we 
ease the IWP design requirement 5 which is concerned with actuator tracking. Also, we 
do not consider the signal conditioning and error weighting blocks since this is just a 
theoretical model and has no practical significance. The CAW gains of this system are 
obtained by first optimizing the J 2 norm of the limited IWP system, and then minimizing 
the structured singular value with a constraint on the J : norm. A 10% deterioration in the 
rms performance was allowed to reduce the directional sensitivity of the system. 

The response of this CAW system to the step input [0.61 0.79] T is presented in 
Figure 2.11. The oscillations exhibited by CAW system are similar to that of the 
unprotected system, but the amplitudes of the oscillations are lesser. The CAW scheme, 
thus maintains the stability and provides a performance better than the unprotected 
system when the system encounters operational limits. This is supported by the structured 
singular value (p) analysis of the limited system. This analysis assumes that the actuator 
limits enter the system as diagonal input uncertainties. The structured singular value is 
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computed by assuming a 10% magnitude bound for the uncertainties in the actuator 
limits. The p plot presented in Figure 2.12 shows that the maximum p of the CAW 
system is less than unity for all diagonal uncertainties satisfying the above magnitude 
bound in the frequency range (0.001 -10)Hz. Hence, the CAW system is robustly stable to 
direction changes of the control vector. 

MAW scheme is implemented for this system in Figure 2.10. The design variable 
P was set at ‘1000’ for this problem. A smaller value for P does not provide adequate 
IWP. It must be noted here that with a high p, the resulting system has a high bandwidth 
and this might cause difficulties in a practical situation. Since this is just an example, we 
assume that it is permissible to have a high P and that there is no limit on the actuator 
bandwidth. The response presented in Figure 2.11 shows that the oscillations and 
overshoot associated with windup are eliminated completely with the MAW scheme. 
Since the MAW scheme maintains the direction of control vector when the system 
encounters operational limits, the performance of the MAW system matches closely with 
the nominal system performance. 

Thus we can conclude that while the limited system A is robustly stable to control 
vector direction changes caused by actuator limits, its performance is sensitive to control 
vector direction changes. Since CAW scheme does not maintain the direction of control 
vector, it causes performance degradation when the actuator limits are encountered The 
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MAW scheme maintains the control vector direction and hence provides a performance 
close to nominal performance when limits are encountered. 


Example 2 : 

The second example is given by : 
P = PoP\ 


with 


4(0.1 + s) | 

Po = — R _1 , P i = 


2(10-5) 
10 + 5 


0 

5-5 ’ 
5 + 5' 


and K = 


1 


4(0.1 + 5 ) 


The response of the nominal and limited systems for a step input of [0.36 0.93 ) 1 
are presented in Figure 2.13. As seen in the figure, the system becomes unstable upon 
encountering actuator limits. It was found that CAW scheme could never stabilize the 
limited system for any set of gains. An unstable response of the CAW scheme is 
presented in Figure 2.14. This is supported by the structured singular value analysis. A 
10% uncertainty in the actuator limits was assumed similar to the previous example. The 
structured singular value plot presented in Figure 2.15 shows that the since the maximum 
p is greater than unity, and hence the CAW system is not robustly stable to uncertain 


actuator limits. 
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The response of the MAW system is presented in Figure 2.13. Since the MAW 
scheme maintains the control direction, it provides a stable system with an acceptable 
performance. We can conclude the following from the above examples : 

• When limited, system A is robustly stable to control vector direction changes. Hence 
the CAW scheme provides a stable solution when actuator limits are encountered. 
However, the response provided by CAW scheme is degraded while MAW scheme 
provides an acceptable performance. This could be due to the sensitivity of the 
performance to the control vector direction. 

• System B when limited, is not robustly stable with control vector direction changes. 
Hence CAW scheme could not stabilize the limited system, while the MAW scheme 
provides a stable, acceptable response. 

• Requirement # 6 presented in section 2.2 has to be satisfied to provide acceptable 
IWP characteristics. While the MAW scheme satisfies this by maintaining the control 
direction, with CAW scheme this is achieved by including the directional sensitivity 


in the performance index definition. 























Chapter 3 


Optimization Technique - GA with Stochastic Coding 


Many of the modem day problems in engineering, medicine, economics, 
sociology, and other areas can be solved in different ways. The question then is to 
determine the strategy that solves the problem in some optimal sense. This is answered by 
quantifying the criteria we are trying to optimize, and then using an optimization 
technique to optimize that criteria. The optimization technique used to solve a non-linear 
optimization problem could play a crucial role in the overall solution process. In some of 
these problems, the optimization step takes the maximum computational time and effort 
expended to solve the problem. Hence, by a judicious choice of the optimization 
technique it might be possible to significantly reduce the overall effort required to arrive 
at an optimal solution. In this chapter, we examine a modified Genetic Algorithm 
optimization technique to obtain the optimal parameters for a test problem given in 
Reference [16]. The modified genetic algorithm approach uses a stochastic parameter 
encoding to find the optimal set of parameters [9], This test problem has non-linear 
interactions between the parameters (i.e. it cannot be solved for each parameter 
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independently), and several locally optimal solutions. The modified GA technique is 
shown to work effectively for this test problem. 


3.1 Conventional Optimization Techniques 

Traditional optimization techniques can be classified as calculus based, 
enumerative, or random search procedures. The calculus based approach could be either a 
direct or an indirect method. In the indirect method, the gradient of the cost function 
(with respect to the parameters being optimized) is set to zero and the resulting equations 
are solved to obtain the optimal set of parameters. In a direct approach, we move in the 
direction of the steepest gradient (of the cost function with respect to the parameters 
being optimized) to find a local optimum. The greatest disadvantage of either of these 
methods is that it is local in scope. It just gives the best solution in the neighborhood of 
the point being considered in the search space. Also, in cases where the gradient 
information is hard (or impossible) to obtain, these methods would fail. 

In an enumerative method, the objective function is evaluated at a large number of 
points in the search space and the optimal solution is the best cost obtained over all the 
points. Such a technique would obviously become inefficient for large search spaces. In a 
random search technique, the objective function is evaluated at some random points in the 
search space. Though this method might provide a global solution, it requires a large 
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number of points to be sampled before the optimal solution is found. Another drawback 
of the random search procedure is that it is difficult to incorporate a stopping criterion in 
the algorithm. 

Hence, we can summarize by stating that the calculus based methods, though 
directed towards an optimum, suffer because of their local scope. The random search 
technique though has the potential to obtain the global optimum, suffers due to the lack of 
a directed search. A search technique referred to as Genetic Algorithms (GA), overcomes 
the shortcomings of the conventional techniques listed above while attempting to retain 
the positive aspects of both of them. In other words, GA is a search technique that 
provides a directed random search [ 17 ]. 


3.2 Genetic Algorithms (GA) 

Genetic Algorithms are search techniques based on principles of natural selection 
and genetics. GAs combine a Darwinian survival-of-the-fittest strategy with a 
probabilistic information exchange between various feasible solutions. GAs are different 
from the conventional gradient based search techniques in the following ways [ 17 ]: 
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1 . GAs operate on a coding of the parameter set, not on the parameters themselves. 

2. GAs search for the optimum solution from a population of points, not from a single 
point. 

3. GAs need information about the objective or cost function, but do not need any 
derivative or other auxiliary information. 

4. GAs use probabilistic transition rules to move across the search space to obtain an 
optimal parameter set. 

Since a GA does not require any derivative information to guide their search, it is 
suitable for problems where the derivative information is difficult or impossible to obtain. 
Other features of the GA that makes it a powerful optimization technique are its 
capability to incorporate apriori knowledge of the solution space, multiple objectives and 
multiple solutions. A typical GA search process is depicted below: 

Step 1 : (Initialization) Obtain a random initial population consisting of n individuals 
representing n points in the search space. Each individual is a binary coded string of the 
parameters of the optimization problem. This binary coded string is referred to as the 
genotype while the actual values of the parameters are referred to as the phenotype. 

Step 2 : (Fitness Evaluation) Each individual of the above population is first decoded to 
obtain the phenotype represented by the individual. These decoded parameter values are 
used to evaluate fitness of each individual using the fitness function. 
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Step 3 : (GA Selection) The selection procedure identifies the fit individuals of the 
current population that would be the parents for the next generation. Two individuals are 
chosen with probabilities proportional to their relative position in the current population, 
measured either by their contribution to the mean fitness value of current generation 
(proportional selection), or by their rank (linear ranking selection). 

Step 4 : (GA Recombination) Two different offspring are produced due to recombination 
of the two parental genotypes of step 3 by means of a crossover with probability p c . Steps 
3 and 4 are repeated until we obtain n individuals for the next generation. 

Step 5 : (Mutation) The offspring obtained in the above step are finally mutated with a 
small probability (p m ). The mutations are assumed to work on individual binary bits, 
either by reversing a one to a zero or otherwise. 

The central processing power of the GA arises from the successful sampling and 
recombination of low order, highly fit, short defining length schemata into strings of 
better fitness. Holland’s schema theorem proves that such low order, highly fit, short 
defining length schemata, also known as building blocks, grow exponentially in a 


population [17J. 
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GAs have been shown to do well on problems with small number of parameters to 
optimize, and are difficult to solve using traditional optimization techniques [9], 
However, for a problem with large number of parameters, large string lengths could result 
when the parameters are coded in a binary form. In reference [9], it has been shown that 
the probability of finding higher order schemata in a population of large strings 
approaches zero rapidly, as the string length increases. This means that the GA would 
attempt to solve the problem by mostly manipulating the lower order schemata, in which 
case, it would take a long time before the optimal solution is found. It is possible to 
increase the probability of the higher order schemata by having a larger population, but 
this calls for enormous computational resources. Hence, we need to modify the GA 
approach so that its effectiveness is maintained while solving problems with large 
parameters. 

Reference [8] examines the techniques to modify the regular GA to maintain its 
effectiveness while solving large parameter problems. In the sensitivity based method, the 
problem of large string lengths is avoided by choosing a few sensitive parameters and 
optimizing the performance index with respect to these sensitive parameters. However, 
this approach is highly problem dependent, and for large parameter problems, even the 
task of determining the sensitivity of each parameter could result in a significant 
computational overhead. Also, if most of the parameters of the optimization problem are 
highly correlated, this approach would fail. 



44 


In GA-local search hybrid technique, the GA searches the solution space at a 
global level while the local/direction search algorithm searches locally around the 
solutions provided by the GA. This enables us to have a coarser encoding structure for the 
GA, thereby reducing the string length. This hybrid technique however results in a 
computational overhead if the local search does not improve the solution provided by 
GA. 


3.3 Stochastic Genetic Algorithms 

The previous sections described the problems faced when we use the GA and the 
modified GA techniques for solving large parameter problems. In this section, we detail a 
GA with stochastic encoding of the parameters (referred to as a Stochastic GA) that 
overcomes these problems and converges quickly to an optimal solution. 

Stochastic Genetic Algorithm is first presented in [9] as an approach to effectively 
solve problems with large number of parameters. Some of the features of the Stochastic 
GA as given in [9] are : 

1 Each discrete possibility as decoded from the binary string, represents a search region 
and not a single value. 

2. The regions defined above are dynamic, (i.e.) the same genotype could represent a 


different region at a different time. 
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3. The search regions are altered based on the GA evolution. 

4. The search region is not explicitly constrained. 

5. No region is completely discarded. 

These features can be implemented in the GA by encoding the search region as a 
binary string. The search region is represented by a multivariate Gaussian distribution 
with a mean vector(p) and the variance matrix(I). The mean vector gives the expected 
values of the parameters in the search region, and the variance matrix gives the 
probability of finding an optimal parameter set in a particular area of the search region. 
The stochastic children are obtained by sampling this multivariate Gaussian distribution. 
In this study, two variations of the stochastic GA were implemented for a test problem. 
The two approaches differ in the way the stochastic children are obtained from the 
parent’s phenotype. In the first approach (Approach A), the stochastic child is obtained 
by varying all the parameters in the parent’s phenotype using the multivariate Gaussian 
distribution (Reference [9]). In this study, the variance matrix used in the multivariate 
Gaussian distribution of Reference [9] is adapted continuously as the GA population 
evolves. This adaptation of the variance matrix helps to exploit the most promising 
regions as the GA explores the search region. In the second approach (Approach B), only 
one parameter is stochastically varied in each string of the GA population. The difference 
between the two approaches can be intuitively pictured as follows. For the stochastic 
child in the approach A, all the features of the child are slightly different from the parent, 
while for the approach B, the stochastic child retains all but one feature of the parent. 


46 


Thus, in approach A, change in the fitness due to the overall change in features is 
detected, where as in approach B; the contribution of each of the parent’s feature to the 
fitness is determined. Hence, with approach B it is possible to fine tune each feature 
independently to obtain an improved fitness. The algorithm that implements the above 
details is presented below. 

Step 1: (Initialization) An initial population of n individuals, characterized by its 
genotype is randomly generated. Each individual’s genotype is a binary string 
representing a search region of the parameters. The search region is encoded in the binary 
string by means of a multivariate Gaussian distribution with mean vector p and a variance 
matrix £. 

Step 2 : (Stochastic Phenotype Variation) 

Approach A : Each of the n individuals produce m offspring, so that a total of mn new 
individuals are available. The search regions represented by these offspring's are obtained 
by displacing the parent’s mean vector(p), with the variance matrix(E), resulting in a new 
mean vector for the offspring. 

Approach B : In this approach, only one parameter in the parent’s mean vector (p) is 
displaced by the variance of that parameter choice. Thus, in both the above approaches, 
the phenotype of the descendent is thus slightly different from that of the parent, while 


the genotype is the same. 
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Step 3 : (Filtering) Out of mn individuals in step 2, only n individuals become parents. 
The phenotypes of the chosen individuals are used to redo the coding, resulting in a 
modified mean vector. In approach A, the variance matrix of the offspring is altered 
based on (1/5)^ success rule. According to this rule, the variance of the Gaussian 
distribution is decreased if atleast one out of five phenotype variations (of the same 
genotype) in step 2 results in an improvement of the performance index. Otherwise, the 
variance is increased. The algorithm that implements this is given in Figure 3.3. The 
variance matrix is not altered in approach B. 

Step 4 : (GA Selection) Two parents are chosen with probabilities proportional to their 
relative position in the current population, measured either by their contribution to the 
mean performance of the current generation (proportional selection) or by their rank 
(linear ranking selection). In a tournament selection procedure, two best individuals from 
a random number of individuals are chosen for the next generation. 

Step 5 : (GA Recombination) The two parents selected in step 4 are recombined with a 
probability of crossover p c , giving rise to two new parental genotypes. Steps 4 and 5 are 
repeated until we have n new individuals representing the next generation. 

Step 6 : (GA Mutation) The new generation obtained above undergoes a mutation 
operation, where the individual bits of each offspring are mutated (reversed from a one to 
a zero, or vice versa) with a small probability (p m ). 
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Initially when the search is started, a symmetric multivariate Gaussian distribution is 
assumed over the search region. The initial Gaussian distribution for a two parameter 
optimization problem is illustrated in Figure 3.1(a). A two dimensional planform view of 
the Gaussian distribution and the parameter choices are shown in Figure 3.2(a). In 
approach A, as the GA population evolves, the distribution shifts towards more promising 
regions and the variance matrix is altered to exploit the promising regions. This is shown 
in Figure 3.1(b) for the two parameter case, where the Gaussian distribution is shown 
after N generations. Figure 3.2(b) shows the planform view of the Gaussian distribution 
after N generations. 


3.4 Application of Stochastic GA for a Test problem 

The stochastic GA methods were applied to a test problem given in (16]. The 
function f (x, ,x 2 ) given below has strong nonlinear interactions between the variables X! 
and x 2 i.e. the optimal value for one parameter cannot be determined independent of the 
other. Such problems are not easily solved by the random hill climbing technique. The 
nonlinear interaction between the variables of the function can further be increased by 
constructing a composite function of more parameters from this primitive function. The 
function F(x h x 2 ....x n ), can in general be a function of any number of variables and is 
constructed from the primitive function of two variables f(Xj,x,+|), as shown. 
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f{x\,x 2 ) = (* 2 | +x 2 2 ) 025 [sin 2 (50(x 2 1 + x 2 2 )° ') + 1.0] 

n - 1 

F(x lr x 2 —Xn) = /(x ll x„) + Y t f(x i ,x i+ 1 ) *,■ €[- 100 , 100 ] 

<=l 

In this study, the function F(X|,x 2 ,....x n ) is optimized for 51 parameters using the 
GA techniques described before. The optimization was carried out with both the regular 
GA and the modified Stochastic GA approaches. For the regular GA approach, a 
population size of 101 was considered. The parameters were assumed to lie in the range 
[-100,100], 6 bits were used to encode each parameter, a probability of crossover of 0.77 
and a probability of mutation of 0.0077 were used. 

For the stochastic GA, 3 offspring were assumed to be produced at the step 2 of 
the algorithm. An initial range of [-100,100] is assumed, and in the approach A, the 
standard deviation in step 3 is decremented by 5%. Also the increment step for the 
standard deviation was set at 2%. (Refer to algorithm in Figure 3.3 for details). The 
population size was set at 5 1 . The probability of crossover and mutation for the stochastic 
GA were set at same levels used for the simple GA, and tournament selection procedure 
was used for both simple GA and stochastic GA. 

For the stochastic GA in approach B, first parameter was fine tuned in the first 
individual of the population, second parameter in the second individual, and so on. Hence 
for this approach, we must have as many individuals in the GA population as the number 
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of parameters. For problems with large number of parameters, we might have to carry out 
a sensitivity analysis, and then use this technique to determine the optimal values for the 
few sensitive parameters. 

3.5 Results of Stochastic GA search 

The results obtained from the application of stochastic GA and the regular GA are 
presented in Figure 3.4. As seen in the figure, the stochastic GA in approach B converges 
quicker than the regular GA and the stochastic GA in approach A. The stochastic GA in 
approach B also converges to a local optimum closer to the global optimum than the 
regular GA and the stochastic GA in approach B. The drawback of approach B, however 
is that we must have a population size equal to the number of parameters of the 
optimization problem. Hence approach B becomes computationally expensive for 
problems with large number of parameters. The stochastic GA approach A when 
implemented for such problems would converge to the optimum quicker than the regular 


GA method. 


Probability 
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Figure 3.1(a) Initial symmetric Gaussian distribution for a two parameter case. 


Probability 
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Figure 3.1(b) Gaussian Distribution for a two parameter case after GA evolution. 


53 


P2 1 

i 



01 1 






010 



© 



001 


o 


©: 


000 



© 




AAA 

AA I 

A 7 A 

All 



TKH3 — DDl m m 


Figure 3.2(a) Planform view of the Initial Gaussian distribution and binary coded 
parameter choices for a two parameter case. 



Figure 3.2(b) Planform view of the Gaussian distribution and the binary coded 
parameter choices after GA evolution (N generations) 
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PI jjarent <— Fitness (Parent) 

PIbestchild <—ritoess(bestchild) 

If ( PI bestchild better than PI jyarent) 

Phenotyp e(Parent) <— Phenoty pe(bestchild) 
end 

If (Stochastic Phenotype Variation in step 2 improves fitness atleast once in 5 variations) 
decrement Varianc ^(Parent Genotype) by 5% 

else 

increment Variance ( Parent Genotype) by 2% 

end 


Figure 3.3 Stochastic GA in approach A. Algorithm that adapts the variance of the 
Multivariate Gaussian distribution that encodes the parameters as binary strings, 
“fitness” returns the fitness of the individual passed to function, Parent and 
bestchild refer to the parent and the best stochastic child obtained at step 2 of the 
algorithm. 






Chapter 4 


Application to STOVL Aircraft Engine Controller Problem 


The advanced maneuvering capabilities of modem aircraft require the propulsion 
system to play a vital part in the aircraft control. Separately designed flight and 
propulsion control systems are however, inadequate for efficient operation of the aircraft 
with reasonable pilot workload [18]. Hence, an Integrated Flight Propulsion Control 
(IFPC) system that accounts for all the subsystem interactions is necessary to obtain an 
optimal system performance with minimal pilot workload. Such a centralized design 
results in a high-order IFPC system that poses implementation problems. The high-order 
centralized controller is therefore partitioned into lower order airframe and propulsion 
controllers with a specified interconnection structure. The partitioning is such that the 
performance and the robustness characteristics of the assembled partitioned controllers 
match that of the centralized controller [18]. Such a partitioning also makes it easier to 
perform independent testing of the airframe and propulsion subcontrollers. Also, the 
partitioning allows the designer to address the airframe and propulsion system 
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non-linearities separately, and to design the integrator windup protection independently 
for each subcontroller [19], 

Preliminary evaluation of the NASA Lewis linear integrated flight/propulsion 
control system has indicated the need for integrator windup protection on the propulsion 
and the airframe subsystems for large command inputs. In this chapter, the CAW and 
MAW techniques will be applied to provide IWP for the linear model of the engine 
subcontroller that forms the part of IFPC system of a STOVL aircraft. 

4.1 The STOVL Aircraft Engine Control System Model 

The aircraft considered for study is a model of a supersonic STOVL aircraft 
powered by a high performance turbofan engine. A schematic diagram of the aircraft with 
various controls is shown in Figure 4.1(a) [20]. The engine control system of the STOVL 
aircraft is equipped with the following controls : 

• ejectors to provide propulsive lift at low speeds and hover. 

• a vectoring ventral nozzle for pitch control and lift augmentation at low speeds. 

• a two-dimensional Convergent-Divergent (2D-CD) vectoring aft nozzle with after 
burner for supersonic flight. 


compressor fan speed. 


The control variables that control the above parameters are respectively : 

• ejector butterfly angle (ETA, degrees) 

• ventral nozzle area (A78, inch 2 ) 

• aft nozzle area (A8, inch 2 ) 

• fuel flow rate to engine (WF36, lbm/hr) 
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Figure 4.1(b) Schematic of the Engine Control System 


The engine state space model can be represented as ( Figure 4.1(b)): 

= A e X e + B e U e 
Y e =C e X e + D e U e 


Where A e ,B e ,C e ,D e are the system matrices of the engine (given in the appendix) 
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X e is the engine state vector, 


WF36* 


FG9 

A8 

Y = 

FGE 

ETA 

' 1 e 

FGV 

A78 


N2 


FG9 - 

aft nozzle thrust (lbf) 

FGE - 

ejector thrust (lbf) 

FGV - 

ventral nozzle thrust (lbf) 

N2 

compressor fan speed (rpm) 


The state space model of the controller can be written as (Figure 4.1(b)): 


* c = A c X c + B c 


e_ zeng 
N2 


y c = 


C c X c + D c 


e_ zeng 
N2 


where A c , B c , C c , D c are the system matrices for the controller (given in the appendix), 
X c is the controller state vector, 


The vector of performance errors : e_ zeng 


FG9_c - FG9 
FGE c- FGE 
FGV c- FGV ; 
N2 c-N2 


ezeng 

N2 


The inputs to the controller are: uc = 
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The outputs of the controller are: 


WF36 

A8 

f = ETA 
A78 


4.2 Selection of limited Actuators 

The engine model considered above has four inputs, four control variables, and 
four actuators. The four actuators control the fuel flow to the engine, ejector butterfly 
valve angle, ventral nozzle area, and the aft nozzle area respectively. It would be ideal if 
we designed the IWP for all possible actuator limit combinations. However, observations 
at NASA Lewis Research Center have shown that in practice that situation never occurs. 
Hence, to simplify the analysis and to study a more realistic situation, we consider a case 
where two of the four actuators are limited. It was expected that the fuel flow to the 
engine would be strongly coupled with the ejector thrust generated by the engine. Since it 
was desired to study the effectiveness of the IWP schemes for systems with strong 
coupling, the fuel flow to engine (WF36) and ejector thrust (FGE) were assumed to be 
limited. This choice of actuator limits is shown in the engine closed loop block diagram 
in Figure 4.2. Here the ejector butterfly valve angle and the fuel flow rate commanded by 
the controller are limited by the hard limits ‘ ujiml ’ and ‘ u_lim2 ’ (blocks 31 and 41). 
Other actuators, the aft-nozzle area and the ventral nozzle area are assumed to be 


unlimited. 
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4.3 IWP Optimization Structure 


For the selected actuators, the IWP performance definition structure was setup as 
shown in the second chapter. The performance definition structure for the engine 
controller problem is given in Figure 4.3. Blocks 1-4, 11-13 are the command loop 
shaping blocks that scale the input magnitude and the frequency spectra. More 
specifically, block 1 1 is the scaling block while the blocks 1-4, 12, 13 are the frequency 
shaping blocks with a first order lag. Blocks 5-10 and 15 weight the performance errors 
and actuator position errors as shown. Block 15 is the scaling factor that is same as the 
inverse of block 1 1, blocks 6-9 weight the performance errors between the IWP system 
and the nominal system, and blocks 5, 10 weight the actuator position errors. The errors 
between the nominal and the IWP systems are obtained in block 14, “Error with IWP”, 
the details of which are shown in Figure 4.4. 


In Figure 4.4, block 7 is the nominal engine control system without limits, while 
block 23 is the limited engine system for which the IWP has to be designed. The 
performance errors and the actuator position errors are shown in this figure. The details of 
block 6 are shown in Figure 4.5, where block 13 is the engine system, block 6 is the 
controller, and block 21 is the scheduling gain to extend the 80 knot nominal design 
controller for the 100 knot engine. The IWP is built inside the control system of the 
limited engine model, i.e. within block 6 of Figure 4.5. 
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4.4 CAW Scheme for Engine Control System 

The CAW scheme is implemented in Figure 4.6. The terms of the gain block 
‘ cengiwp ’ are the gains of the CAW scheme and are to be obtained by optimizing the 
performance index. The actuator error is fed to the ‘ ceng iwp ' block as shown and the 
outputs from the ‘ ceng iwp ' block are added to the state derivative calculations. With 
CAW implemented, equations for the controller can be modified as : 

e_ zeng 

X c = A c X c + B c N2 

e_ ueng 

e_ zeng 

Y c = C c X c + D c N2 

e_ ueng 

e_ETA 

where e ueng = is the vector of actuator errors. 

e_Wr36 

With zero actuator error, the IWP gains do not contribute to the state derivative 
calculations and hence the CAW scheme is memoryless. 

4.5 MAW Scheme for Engine Control System Model 

The implementation of the MAW scheme in the engine controller is shown in 
Figure 4.7, where the blocks 3 and 6 implement the MAW scheme in the nominal control 
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system as shown. During the optimization stage, a value for a (windup factor) is assumed 
and the design variable p is obtained for the best system performance. With the MAW 
scheme implemented, the controller equations can be given as: 


X , 


4 


) X +B 

ezeng 

) c c 

N2 


ezeng 
N2 

where e zeng is the vector of performance errors. When the actuators are not limited, 
a=l, and we get back the nominal control system. Thus, the MAW scheme is also 
memoryless. The system matrices of the controller, [A c , B c , C c D c ] that appear in Figure 
4.7 are listed in the appendix. 


T c = C c * c+ aD c 


4,6 Optimization of the Performance Index 


The algorithm to optimize the 2-norm (rms norm) error of the IWP optimization 
structure is presented in Figure 4.8. First, we ensure that the chosen set of gains stabilize 
the controller and the closed loop IWP system. If the IWP gains form a stable system, the 
2-norm of the system is obtained from “rms” function available in MatrixX software. The 
IWP gains that yield the minimum 2-norm are obtained using the stochastic GA 
optimization technique described in chapter 3. 
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The algorithm to optimize the performance index of the CAW system is given in 
Figure 4.9, where the structured singular value (p) of the IWP system is minimized with a 
constraint on the 2-norm. The structured singular value (ssv) of the IWP system is 
obtained from “ ssv ” function available in the MatrixX software. A 10% bound for the 
diagonal uncertainties and a frequency range (0.001-10) Hz are assumed for the p 
analysis. ‘ k ’ is the factor (>1) that determines the 2-norm performance degradation we are 
willing to tolerate in order to improve the structured singular value of the system. In this 
study, ‘k’ is assumed to be 1.1, which means that we allow a 10% degradation in the 2- 
norm performance to minimize directional sensitivity of the IWP system. First few steps 
of this algorithm ensure that the IWP gains form a stable controller and a stable closed 
loop IWP system. If a set of gains do not yield a stable solution, or if they do not satisfy 
the 2-norm constraint, a high penalty is imposed on them. Thus, the parameters that 
satisfy the rms norm constraint and minimize the structured singular values of the system 
are obtained as the gains for the CAW scheme. In this algorithm, J 2mm is the minimum 2- 
norm of the IWP optimization structure given in Figure 4.3. 

The algorithm to compute the performance index for the MAW scheme is similar 
to that shown in Figure 4.8. In the MAW approach, we need not include the structured 
singular value computations. This is because, the MAW method preserves the direction of 
the control vector and hence there is no directional uncertainty. Thus the performance 
measure for the MAW scheme is just the 2-norm of the optimization structure presented 
in Figure 4.3. Also in the MAW approach, we have just one design variable ((3) to be 
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determined for optimum performance. The rms performance measure was found to 
monotonically decrease with p, indicating that a high value of p is ideal for good 
performance. However, with high beta, the system bandwidth increases and also, it may 
not be practically feasible to implement a high p. Hence, from these considerations the 
value for P was chosen to be 25 for this study. 


4.7 IWP implementation Results for the Engine Controller Problem 

In this section, we present the results of the application of IWP techniques for the 
engine controller of a STOVL aircraft and compare the system performance obtained 
using the CAW and the MAW schemes. 

4.7.1 Nominal and limited system response 

During simulation of the engine control system, a FGE step input command of 
1000 lbf and a N2 step input command of 100 rpm were assumed. The nominal system is 
the linear control system without any actuator limits. For the limited control system, fuel 
flow to the engine and the ejector butterfly valve angle were assumed to be limited. The 


limits on ETA and WF36 were set at 7.0 degrees and 450 Ibm/hr respectively. 
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The FGE command tracking response of the nominal and the limited system 
(without IWP) is presented in Figure 4. 10. The effects of integrator windup of the limited 
system are clearly seen in this figure. The limited system exhibits oscillations, high 
steady state error, and a degraded performance when the step command is terminated 
after 5 seconds. The degraded performance following the termination of the step 
command is because, the actuators of the limited system do not get off the limits until the 
integrators unwind. The ejector angle commanded by the controller to track the FGE step 
command is given in Figure 4.1 1. As seen from this figure, the nominal control system 
requires an ejector angle of about 9 degrees to track the FGE step command. When the 
ejector angle is limited at 7 degrees, the controller’s integrators windup resulting in poor 
tracking of the limited actuator. 

The N2 command tracking response of the nominal and the limited system 
(without IWP) is presented in Figure 4.12. Again, the limited system exhibits increased 
oscillations and sluggish behavior (increased rise and settling times). The fuel flow 
commanded by the nominal and the limited controllers to track the N2 command is 
presented in Figure 4.13. The nominal system requires about 570 lbm/hr of fuel to track 
the N2 command. When the fuel flow is limited at 450 lbm/hr, the windup causes poor 


actuator tracking as shown. 
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4.7.2 CAW scheme for Engine Control System - Results 

The CAW gains were obtained by first minimizing the 2-norm of the optimization 
structure (J 2m j n ), and then minimizing the structured singular value of the IWP system 
with a constraint on the 2-norm in terms of ^min- Structured singular value analysis for 
the CAW gains is presented in Figure 4.14. Since the maximum p obtained in Figure 4.14 
is less them unity, the CAW system is robustly stable for changes in control vector 
direction. 

When the CAW scheme is implemented to provide IWP, the FGE command 
tracking response of the CAW system matches closely with that of the nominal system as 
shown in Figure 4.10. The CAW scheme eliminates the oscillations and the high steady 
state error exhibited by the unprotected system. When the step command terminates, the 
system with IWP gets off the limits quickly, and tracks the FGE almost exactly as the 
nominal system. Also, when IWP is implemented, an accurate tracking of the limited 
actuator is achieved. This is shown in Figure 4.11, where the ETA commanded by the 
controller is at the limit (7 degrees). 

The command tracking response of the CAW system for the N2 command is 
shown in Figure 4.12. The CAW system exhibits reduced oscillations, a smaller steady 
state error and is less sluggish than the unprotected system. The IWP controller tracks 
limited WF36 accurately as seen in Figure 4.13. We can also observe that the IWP 
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controller demands a higher fuel flow than the nominal system during the time period 
(2-5) seconds. This is due to the limit on the ejector angle, ETA. The control variables of 
this system are all coupled and hence, when one actuator is limited, the controller 
attempts to track the commands with an increased value for other variables. Hence, when 
the ejector angle is limited, the control system attempts to track the FGE command with 
an increased fuel flow to the engine. 


4.7.3 MAW scheme for Engine Control System - Results 

The results of the implementation of the MAW scheme for the engine controller 
are presented in Figures 4.10 and 4.12. Figure 4.10 shows the FGE command tracking 
response of the MAW system. Figure 4.12 shows the N2 command tracking response of 
the IWP system with MAW scheme. We can observe from these figures that with MAW 
scheme, the system gets off the limits quickly and provides a response close to the 
nominal system when the step command is terminated. However, MAW scheme 
increases the rise time and results in a high steady state error for this system as seen in 
Figures 4.10 and 4.12. Figures 4.11 and 4.13 indicate that with MAW scheme, an 
accurate actuator tracking is attained when the system encounters the operational limits. 
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4.7.4 Comparison between CAW and MAW schemes 

In Figures 4.15 and 4.16, the steady state error variations on FGE and N2 
commands with ETA and WF36 limits are presented. Figure 4.15 shows the variation of 
the steady state error on the commanded FGE as the limit on ejector angle is increased 
from 1 degree to 6 degrees. For this case, input to the system is a FGE step command of 
1000 lbf, and a limit of 450 lbm/hr is set for fuel flow. The CAW system results in a high 
steady state error when stringent limits are imposed on ETA. As the limit on ETA is 
relaxed, the error reduces initially but stays constant for ETA greater than 4 degrees. The 
error does not decrease further due to the limit on fuel flow. For a given fuel flow rate to 
the engine, there is a maximum ejector thrust that could be generated. As the ejector 
valve angle is increased from zero, the ejector thrust increases until the maximum value is 
attained. Once the maximum is attained, FGE can be increased further only by increasing 
the fuel flow to the engine. Without the limits on fuel flow, the control system tends to 
track the FGE command accurately for any limited ETA (>0) by appropriately increasing 
the fuel flow to the engine. Figure 4. 1 5 shows that the steady state error is higher with the 
MAW scheme for any ETA limit. 

Figure 4.16 shows the variation of the steady state error on N2 command as the 
limit on WF36 is varied between 200 lbm/hr to 600 lbm/hr. For this case, a N2 step 
command of 100 rpm is the input and ETA limit is set at 7.0 degrees. As seen in the 
figure, MAW system results in a higher steady state error than the CAW system. With 
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CAW scheme, steady state error approaches zero as the limit on fuel flow is relaxed. 
However, with MAW scheme the errors attain a constant value and do not approach zero 
even if the limits on fuel flow are removed. We can conclude from these results that for 
the engine controller problem, with the fuel flow and the ejector valve being limited, the 
CAW scheme provides better anti-windup protection than the MAW scheme. 



Rudder 




nozzle 


Figure 4.1 Control Actuators of a STOVL aircraft. 
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Figure 4.2 Limited Engine System. 





Figure 4.3 IWP Optimization Structure. 











Nominal System 
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Figure 4.4 Error with 1WP System. 
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Figure 4.5 IWP Closed Loop System 
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igure 4.6 CAW Scheme Controller. 
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Figure 4.7 MAW Scheme Controller. 
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While (Solution not converged) 
cawjwp 4- IWP gains 

If ( engine controller and closed loop system ) are stable 
J 2 <— rms(“IWP Optimization Structure”) 

else 

J 2 <— high penalty 
end 

endwhile 

Figure 4.8 Algorithm for RMS norm optimization, “rms” is the function that 
returns the 2-norm of the system. 

J 2mm minimum rms. norm of “IWP optimization structure”. 

While (Solution not converged) 

cawjwp <— IWP gains chosen 

If ( engine controller and closed loop system) are stable 
J 2 <— rms(“IWP Optimization Structure”) 
if (J 2 < k*J 2 mm) 

p 4— ssv(“IWP Engine System”) 

else 

p 4— high penalty 

end 

else 

p 4- high penalty 
end 

endwhile 

Figure 4.9 Algorithm for CAW gain optimization, “rms” returns the 2-norm of the 
system passed to the function, and “ssv” returns the maximum structured singular 
value of the system passed to the function. 


ETA (Degrees) n FGE Thrust (Ibf) 
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Nominal System 

Limited System without I WP 

— . — CAW Scheme 

— — MAW Scheme 



4.10 FGE command tracking response for nominal, limited (without IWP), 
CAW and MAW Control Systems. 



Figure 4.1 1 Ejector Angle (ETA) commanded by the nominal, limited, CAW and 
MAW controllers. The limit on ETA is 7 degrees. 




80 


Nominal Sytstem 
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— • — CAW Scheme 

— — MAW Scheme 
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Figure 4.12 Nominal, limited system (w ithout IWP), CAW and MAW Control 
System responses for a N2 step command of 100 rpm. 
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Figure 4.13 Fuel flow (WF36) commanded by the nominal, limited, CAW and 
MAW controllers. The limit on WF36 is 450 Ibm/hr. 
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ETA Limit (deg) 

Figure 4.15 Steady state error on FGE step command of 1000 Ibf with varying 
ETA limits. Fuel flow (WF36) is limited at 450 Ibm/hr. 
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Figure 4.16 Steady state error on N2 step command of 100 rpm with varying limits 


on fuel flow. ETA limit is set at 7.0 degrees. 



Chapter 5 


Conclusions and Recommendations 

Two anti-windup compensation techniques based on the back-calculation 
principle, namely the Conventional Anti-Windup (CAW) and the Modified Anti-Windup 
(MAW) techniques were described and implemented for three problems. An unified 
performance index has been defined for the CAW scheme such that the directional 
sensitivity of the CAW system is minimized while providing an acceptable performance. 
Since the MAW scheme maintains the control direction, we just consider the 2-norm 
performance measure to obtain an optimal performance for the MAW system. The results 
obtained from the IWP implementations for these problems are summarized below. 

In the first illustrative example (System A) considered in chapter 2, the limited 
system was robustly stable to changes in the control direction. Hence the CAW scheme 
could provide a stable IWP controller, but the performance of the CAW system was 
deteriorated. This could be due to the sensitivity of system performance to the control 
direction. Since the MAW approach maintains the direction of the control vector, a stable 
system with an acceptable performance was obtained with the MAW system. 
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In the second example (System B) of chapter 2, the limited system was not 
robustly stable for changes in control vector direction caused due to actuator limits. 
Hence it was not possible to obtain a CAW controller that stabilized the limited system. 
The MAW scheme provided a stable controller with an acceptable performance. 

The limited engine controller system considered in chapter 4 was robustly stable 
for uncertain actuator limits, and the simulation results indicate that CAW system 
provides a performance better than the MAW scheme. Since in the MAW scheme, the 
design variable p is the only choice to obtain an optimal performance, it is not always 
possible to tune the system to a desirable performance. Based on these results, we can 
draw the following conclusions : 

• The directional sensitivity of the IWP system must be included in the performance 
index to minimize the system sensitivity to control direction changes caused when 
actuators encounter operational limits. 

• If the limited system is robustly stable to diagonal uncertainties, then CAW approach 
would provide a stable IWP system. The CAW scheme would perform better than the 
MAW scheme if the performance of the system is not very sensitive to the control 
direction. However, if the performance is sensitive to the control direction, then 
MAW scheme would perform better. 
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If the CAW system is not robustly stable to control direction changes, the MAW scheme 
could provide adequate windup protection. 

The following recommendations are made based on the results obtained in this 
study. Given a linear control system for which a IWP technique is to be implemented to 
satisfy the general IWP requirements, we first attempt to implement the CAW technique 
for the system. If it is possible to obtain a stable closed loop IWP system with the CAW 
scheme, we optimize the 2-norm performance of the CAW system after including 
appropriate loop shaping and frequency weighting blocks in the performance definition 
structure. Once the set of CAW gains that minimize the 2-norm of the IWP system is 
obtained, we can minimize its directional sensitivity by minimizing the structured 
singular value of the IWP system with a constraint on the 2-norm. The CAW gains thus 
obtained will minimize the directional sensitivity of the closed loop system while 
ensuring an acceptable performance. 

If a set of CAW gains that stabilize the limited system cannot be obtained, the 
limited system is not robustly stable to control vector direction changes. Hence we must 
resort to the MAW scheme to maintain the control direction when the actuator encounters 
the operational limits. 

As a general conclusion, if the CAW scheme optimized with the unified 
performance measure shows robust stability, it is preferred over the MAW scheme. 
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Appendix 


The system matrices for the Engine are given by : 



- 4.7000 10e + 00 

3.018767e + 00 

4.967368e-0l 

1.602869e-01 .. 


4.065682e-01 

-3.244523e + 00 

6.536386e-0l 

2.265022e-01 . 
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- 1 958345e-01 
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7.904529e- 0 1 
1 278608e-03 
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I .OOOOOOe + 00 
1. 64265 le + 00 
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- 1 425019e- 17 
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-3.4l7969e-02 
- 1.379395e-02 


-1.122547e-22 ... 
-2.856445e-02 ... 
-2.563477e-02 ... 
-9.948730e-03 ... 


-8.348264e-24 

1.5869l4e-02 

1.757813e-02 

6.225586e-03 


1 .3857 1 9e - 23 
8.91 1 1 33e - 02 
8.666992e-02 
3. 1921 39e-02 


Dr 


-8.938996e- 19 
l . 1 99692e - 0 1 
1 . 1 40380e - 0 1 
4.239934e-02 


3.256075e-20 
9.158529e + 00 
- 1.245736e + 0l 
-4.720947e + 00 


1.049588e - 18 
-2.332959e + 01 
8.738379e + 0l 
-8.239990e +00 


7.984864e-2l 

- 1.277783e + 01 

- 1.189714e + 01 
1 .392403e + 0 1 


O.OOOOOOe +00 
-2.36621 le + 01 
-2.202637e + 01 
-8.377686e + 00 


The system matrices for the Engine Controller are given by : 


A c 


-1.086374e-02 
4.886350e-05 
- 1.365098e-05 
O.OOOOOOe +00 
O.OOOOOOe +00 
O.OOOOOOe +00 
O.OOOOOOe + 00 


4.l05142e-05 
- 1 .073942e -02 
1.021952e-05 
O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe +00 


-6.839382e-06 
7.6897 lOe- 06 
- l.071602e-02 
O.OOOOOOe +00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 


O.OOOOOOe + 00 ... 
O.OOOOOOe + 00 ... 
O.OOOOOOe + 00 ... 
-7.25l576e-03 ... 
-3.67l925e-03 ... 
- 6.45342 le- 03 ... 
6.830768e-03 ... 


O.OOOOOOe + 00 O.OOOOOOe +00 O.OOOOOOe + 00 
O.OOOOOOe + 00 O.OOOOOOe + 00 O.OOOOOOe + 00 

O.OOOOOOe + 00 O.OOOOOOe + 00 O.OOOOOOe + 00 
3.4 1 3 1 43e - 04 1.879533e-03 3.298592e-04 

-1.451 197e-02 -3.628540e-02 -6.2371 1 8e - 02 
-3.853371e-02 - 1.295285e-01 -3. 1 14797e-01 
5.258483e-02 2.589395e-01 -2.926340e + 00 
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Be 


4.93023 le-04 
-2.493989e-03 
4.544243e-04 
O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 


7.559844e-04 
3.3705 12e -04 
9.250075e-04 
O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 


2. 145 1 38e - 03 
8.464059e-05 
-1.405475e-03 
O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 


O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 
- 1.0l4878e-02 
-2.561923e-03 
-4.5013 12e-03 
4.778705e-03 


O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 
4.633683e-06 
-2.08l017e-04 
- 7.6027 18e- 04 
-2.599835e-0 3 


C c 


1.91 2883e + 03 
-2.96541 le + 01 
-2.33 1950e + 00 
-4.905758e + 01 


-6 844580e + 02 
-3.75 1226e + 01 
1 01 8938e + 01 
4.66 1629e + 01 


1 .33677 1 e + 02 
l.278786e+0I 
1.56445 le +01 
-5.2358 1 le + 01 


-1.668799e + 03 ... 
-6.990420e + 00 ... 
- 5.892 188e + 00 ... 
-3.684842e + 01 ... 


1.606993e + 02 
6.471627e + 00 
1 .673976e + 00 
- 1 .35039 1 e + 0 1 


4.951522e + 02 
-9.979779e + 00 
2.7 1942 le + 00 
-2.139078e + 01 


2.779248e + 02 
-8.76634 le + OO 
4.858876e + 00 
-3.642707e + 00 


Dc 


O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe +00 


O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 


O.OOOOOOe + 00 
O.OOOOOOe +00 
O.OOOOOOe + 00 
O.OOOOOOe + 00 


1.242878e + 00 
3.504186e-05 
- 9.8073 13e- 04 
-4.619351e-03 


-3.416672e + 00 
- 3.268677e-02 
-2.097606e.02 
-9.870098e-02 


The scheduling gains for the 100 knot engine model to extend the 80 knot nominal design 
controller is given by the matrix : 


9. 1 1 8556e - 01 
-4.9997 17e- 01 
-1.461968e - 01 
1 740255e- 01 


- 1 847999e - 01 
2.278030e + 00 
2. 1 13056e-01 
- 1. 604945c- 01 


6. 108760e - 02 !.466123e-01 

4.999329e - 01 4.999664e-0l 

1.250743e + 00 - 1.722879e-03 

- 1 629085e - 0 1 6.807980e - 0 1 



